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1. INTRODUCTION 

The idea of neutrino oscillations was discovered by Bruno Pontecorvo in the late 50s in 
analogy with K^-K^ oscillations [1, 2]. In essence, neutrino oscillations are lepton flavor 
transitions which depend on the distance and time of neutrino propagation between a 
source and a detector. This is a quantum-mechanical effect due to neutrino mixing, 
i.e. the fact that flavor neutrinos are coherent superpositions of massive neutrinos. The 
oscillations are caused by the interference of the different massive neutrinos, which have 
different phase velocities. 

Since in the late 1950s only one active flavor neutrino was known, the electron 
neutrino, Pontecorvo invented the concept of a sterile neutrino v^- [3], which does not 
take part in weak interactions. The muon neutrino was discovered at Brookhaven in 1962 
in the first accelerator neutrino experiment of Lederman, Schwartz, Steinberger, et al. 
[4], following the independent feasibility estimates of Pontecorvo [5] and Schwartz [6]. 
Since then, it became clear that oscillations between different active neutrino flavors are 
possible if neutrinos are massive and mixed^ Indeed, in 1967 Pontecorvo [3] discussed 
the possibility of a depletion of the solar Vg flux due to Vg — > (or Vg v^) transitions 
before the first measurement in the Homestake experiment [9]. In 1969 Gribov and 
Pontecorvo [10] discussed solar neutrino oscillations due to Vg-v^ mixing. 



In 1962 Maki, Nakagawa, and Sakata [7] considered for the first time a model with Vg-V^ mixing of 
different neutrino flavors. Unfortunately, this model did not have any impact on neutrino mixing research, 
since its existence was unknown to the community until the late 70s [8]. 



The standard theory of neutrino oscillations was developed in 1975-76 by Eliezer and 
Swift [11], Fritzsch and Minkowski [12], Bilenky and Pontecorvo [13, 14]. In this theory, 
massive neutrinos are treated as plane waves, having definite energy and momentum. 
Such a description, however, is not completely consistent, because energy-momentum 
conservation implies that the creation and detection of massive neutrinos with definite 
energies and momenta is possible only if all the particles involved in the production 
and detection processes have definite energies and momenta. The problem is that in this 
case energy-momentum conservation cannot hold simultaneously for different massive 
neutrinos and the production and detection of a superposition of different massive 
neutrinos are forbidden. In order to overcome this problem, it is necessary to treat 
neutrinos and the other particles participating in the production and detection processes 
as wave packets, as discussed in section 5. 

The plan of this paper is as follows: in section 2 we review the standard theory of 
neutrino oscillations, highlighting the main assumptions, which are discussed in the 
following sections; in section 3 we discuss the definition of flavor neutrino states; in 
section 4 we present a covariant plane-wave theory of neutrino oscillations; in section 5 
we discuss the necessity of a wave-packet treatment of neutrino oscillations, in section 6 
we present a simple quantum-mechanical wave-packet model of neutrino oscillations, 
and finally in section 7 we draw our conclusions. 

2. STANDARD THEORY OF NEUTRINO OSCILLATIONS 

Neutrino oscillations are a consequence of neutrino mixing: 

yadx) =Y,UakVkL{x) (a = e,jU,T), (1) 
k 

where VaL{x) are the left-handed flavor neutrino fields, Vki{x) are the left-handed mas- 
sive neutrino fields and U is the unitary mixing matrix (see Refs. [8, 15, 16, 17, 18]). 

Since a flavor neutrino Va is created by V^^(x) in a charged-current weak interaction 
process, in the standard plane-wave theory of neutrino oscillations [11, 12, 13, 14, 8], it 
is assumed that is described by the standard flavor state 

\ya) = Y,Kk\yk), (2) 

k 

which has the same mixing as the field V^^(x) . 

Since the massive neutrino states \Vk) have definite mass m;^ definite energy Ej^, 
they evolve in time as plane waves: 

ii-Wkit)) = Ho \Vk{t)) = Ek\Vk{t)) , |v,(0) = e-'^^' |v,) , (3) 



where Hq is the free Hamiltonian operator, 
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Ek ^Pk + f^k, 



(4) 



and \Vk{t = 0)) = \ Vk) (all the massive neutrinos start with the same arbitrary phase). 
The resulting time evolution of the flavor neutrino state in Eq. (2) is given by 

k fi=e,ix,x \ k / 

Hence, if the mixing matrix U is different from unity {i.e. if there is neutrino mixing), 
the state | Va(?)), which has pure flavor a at the initial time f = 0, evolves in time into a 
superposition of different flavors. The quantity in parentheses in Eq. (5) is the amplitude 
of Va transitions at the time t after Va production. The probability of Va — Vjg 

transitions at the time t = T of neutrino detection is given by 



Pa^{T) = \{v^\Va{T))\' = 



pk 



i:U*,,UpkUajU*pje-^^^^-^^^r (6) 

kj 



One can see that Paji{T) depends on the energy differences — Ej. In the standard 
theory of neutrino oscillations it is assumed that all massive neutrinos have the same 
momentum p. Since detectable neutrinos are ultrarelativistic^, we have 



m? Am?, 



E + '-^, Ek-Ej^ — ^, (7) 
2E' " J 2E ' 

where Am|y = ml — mj and E = \p\ is the energy of a massless neutrino (or, in other 
words, the neutrino energy in the massless approximation). In most neutrino oscilla- 
tion experiments the time T between production and detection is not measured, but the 
source-detector distance L is known. In this case, in order to apply the oscillation prob- 
ability to the data analysis it is necessary to express T as a function of L. Considering 
ultrarelativistic neutrinos, we have T L, leading to the standard formula for the oscil- 
lation probability: 



kJ 



^ / Am?,L\ 

P«^(L,£) =£c/*,C/^,[/a,[/^.exp -i-^ ] . (8) 



Summarizing, there are three main assumptions in the standard theory of neutrino 

oscillations: 

(Al) Neutrinos produced or detected in charged-current weak interaction processes are 
described by the flavor states in Eq. (2). 

(A2) The massive neutrino states \Vk} in Eq. (2) have the same momentum {^'equal- 
momentum assumption"). 



^ It is known that neutrino masses are smaller than about one eV (see the reviews in Refs. [19, 20]). Since 
only neutrinos with energy larger than about 100 keV can be detected (see the discussion in Ref. [21]), in 
oscillation experiments neutrinos are always ultrarelativistic. 



(A3) The propagation time is equal to the distance L traveled by the neutrino between 
production and detection (''time = distance assumption"). 

In the following we will show that the assumptions (Al) and (A3) correspond to ap- 
proximations which are appropriate in the analysis of current neutrino oscillation exper- 
iments (section 3 and 5, respectively). Instead, the equal-momentum assumption (A2) is 
not physically justified [22, 23, 24, 25, 26], as one can easily understand from the appli- 
cation of energy-momentum conservation to the production process^. However, in sec- 
tion 4 we will show that the assumption (A2) is actually not necessary for the derivation 
of the oscillation probability if both the evolutions in space and in time of the neutrino 
state are taken into account. 

3. FLAVOR NEUTRINO STATES 

The state of a flavor neutrino Va is defined as the state which describes a neutrino 
produced in a charged-current weak interaction process together with a charged lepton 
£j or from a charged lepton £~ (£^ = e^,/i^,T^ for a — e,ix,T, respectively), or the 
state which describes a neutrino detected in a charged-current weak interaction process 
with a charged lepton £~ in the final state. In fact, the neutrino flavor can only be 
measured through the identification of the charged lepton associated with the neutrino 
in a charged-current weak interaction process. 

Let us first consider a neutrino produced in the generic decay process 

Pl^PF + 4 + Va, (9) 

where Pi is the decaying particle and Pp represents any number of final particles. For 
example: in the pion decay process ;r+ — > ju"*" -I- we have Pi = 7^*", Pp is absent 
and a = /i; in a nuclear jS"*" decay process N(A,Z) — > N(A,Z — \)+e'^ + Vg we have 
Pi = N(A,Z), Pp = N(A,Z— 1) and a = e. The following method can easily be modified 
in the case of a produced in the generic scattering process -|- Pi — > Pp -|- by 
replacing the £^ in the final state with a l'^ in the initial state. 

The final state resulting from the decay of the initial particle Pi is given by 

l/) = S|Pi), (10) 

where S is the 5-matrix operator. Since the final state \ f) contains all the decay channels 
of Pi, it can be written as 

l/)=E-<fc|Vfc,4,PF) + ..., (11) 
k 

where we have singled out the decay channel in Eq. (9) and we have taken into account 
that the flavor neutrino Va is a coherent superposition of massive neutrinos V^. Since 



A different opinion, in favor of the equal-momentum assumption, has been recently expressed in 
Ref. [27]. On the other hand, other authors [28, 29, 30] advocated an equal-energy assumption, which 
we consider as unphysical as the equal-momentum assumption. 



the states of the other decay channels represented by dots in Eq. (11) are orthogonal to 
\Vk,£'^,P-p) and the states |v^,£^,Pf) with different ks are orthonormal, the coefficients 
^^j^ are the amplitudes of production of the corresponding state in the decay channel in 
Eq. (9): 

^ak = (Vfe,4,PF|/) = {ykA,PF\ S |Pi) . (12) 

Projecting the final state in Eq. (11) over |£^,Pf) and normalizing, we obtain the flavor 
neutrino state [31, 32, 17, 33] 




-1/2 

KfelVfc). (13) 



Therefore, a flavor neutrino state is a coherent superposition of massive neutrino states 
|Vyt) and the coefficient .af^j^ of the kth massive neutrino component is given by the 
amplitude of production of v^. Since, in general, the amplitudes depend on the 
production process, a flavor neutrino state depends on the production process. In the 
following, we will call a flavor neutrino state of the type in Eq. (13) a "production flavor 
neutrino state". 

Let us now consider the detection of a flavor neutrino through the generic charged- 
current weak interaction process 

Va+Di^DF + £-, (14) 

where Dj is the target particle and Dp represents one or more final particles. In general, 
since the incoming neutrino state in the detection process is a superposition of massive 
neutrino states, it may not have a definite flavor. Therefore, we must consider the generic 
process 

v + Di, (15) 

with a generic incoming neutrino state |v). In this case, the final state of the scattering 
process is given by 

|/) = S|v,Di), (16) 
This final state contains all the possible scattering channels: 

I/) = |Df, £„) + ..., (17) 

where we have singled out the scattering channel in Eq. (14). We want to find the 
component 

|Va,Di)=£<l^JVfc,Di) (18) 

k 

of the initial state |v,Di) which corresponds to the flavor a, i.e. the component which 
generates only the scattering channel in Eq. (14). This means that |Df,^^) = S | Va,Di). 
Using the unitarity of the mixing matrix, we obtain 

|Va,Di) = S^|DF,£-). (19) 



From Eqs. (18) and (19), the coefficients are the complex conjugate of the amplitude 
of detection of in the detection process in Eq. (14): 

<° = (V;t,Di|S^|DF,4). (20) 

Projecting | v^, Di) over |Di) and normalizing, we finally obtain the flavor neutrino state 
in the detection process in Eq. (14): 

iv°)=(|EKi^j i:<iv.). (21) 

In the following, we will call a flavor neutrino state of this type a "detection flavor 

neutrino state". 

Although the expressions in Eqs. (13) and (21) for the production and detection 
flavor neutrino states have the same structure, these states have different meanings. 
A production flavor neutrino state describes the neutrino created in a charged-current 
interaction process, which propagates out of a source. Hence, it describes the initial 
state of a propagating neutrino. A detection flavor neutrino state does not describe a 
propagating neutrino. It describes the component of the state of a propagating neutrino 
which can generate a charged lepton with appropriate flavor through a charged-current 
weak interaction with an appropriate target particle. In other words, the scalar product 

Aa = (vD|v) (22) 

is the probability amplitude to find a by observing the scattering channel in Eq. (14) 
with the scattering process in Eq. (15). 

In order to understand the connection of the production and detection flavor neutrino 
states with the standard flavor neutrino states in Eq. (2), it is useful to express the 5- 
matrix operator as 

S=l-iJ d^^HccW , HccW = ^jlix)jP{x) , (23) 

where Gp is the Fermi constant (we considered only the first order perturbative contri- 
bution of the effective low-energy charged-current weak interaction Hamiltonian). The 
weak charged current jP (x) is given by 

f{x)= I V^{x)f(l-^)ea{x)+hP{x) 
a=e,n,z 

= E LKkM^)f(l-f)ia{x)+hP{x), (24) 
a=e,iJ.,z k 

where hP (x) is the hadronic weak charged current. The production and detection ampli- 
tudes ^^i. and =2/5 can be written as 



-^ak -^ak^ak^ ^ak-^ak^ak^ 



(25) 



with the interaction matrix elements 



= S / "^'^ ^"'^ ^« I ^^^^ f{i-"^)W |0) 4'^''' (x) , (26) 



V2 

^ak = i^J d\(v,|nW/ (l-r^)^aWK«)/?-^°^*(x). (27) 

Here 7p'^^^(jc) and /p'^^^(x) are, respectively, the matrix elements of the Pj — > Pp and 
Di Dp transitions. 

Using Eq. (25), the production and detection flavor neutrino states can be written as 
|v^)=I^=^— f/^,|v,), (28) 

\^a) = L , Kk I V.) . (29) 



These states have a structure which is similar to the standard flavor states in Eq. (2), with 
the relative contribution of the massive neutrino Vk proportional to The additional 
factors are due to the dependence of the production and detection processes on the 
neutrino masses. 

In experiments which are not sensitive to the dependence of ^^j^ and on the 
difference of the neutrino masses, it is possible to approximate 



In this case, since 

E|f/afe|'=l, (31) 

k 

we obtain, up to an irrelevant phase, the standard flavor neutrino states in Eq. (2), which 
do not depend on the production or detection process. Hence, the standard flavor neu- 
trino states are approximations of the production and detection flavor neutrino states in 
experiments which are not sensitive to the dependence of the neutrino interaction rate 
on the difference of the neutrino masses. All neutrino oscillation experiments have this 
characteristic: since the detectable neutrinos are ultrarelativistic, neutrino oscillation ex- 
periments are insensitive to any effect of neutrino masses in the production and detection 
processes. Therefore, the assumption (Al) in the standard theory of neutrino oscillations 
is correct as an appropriate approximation in the analysis of neutrino oscillation experi- 
ments. 



4. COVARIANT PLANE-WAVE THEORY OF OSCILLATIONS 

In this section we show that the equal-momentum assumption (A2) can be avoided by 
considering not only the time evolution of the neutrino states, as in the standard theory, 
but also their space dependence. 



Let us consider a neutrino oscillation experiment in which Va — ^ V^g transitions are 
studied with a production process of the type in Eq. (9) and a detection process of 
the type in Eq. (14). In this case, the produced flavor neutrino Va is described by 
the production flavor state |v^) in Eq. (13). If the neutrino production and detection 
processes are separated by a space-time interval {T,L), the neutrino propagates freely 
between production and detection, evolving into the state 

|v(r,L)) = ^-E^+'P^|vP), (32) 

where E = Hq and P are, respectively, the energy and momentum operators. This is 
the incoming neutrino state in the detection process. The amplitude of the measurable 
Vce — > transitions is given by the scalar product in Eq. (22): 

Aap{T,L) = {vf\v{T,L)) = {vf\e-'^^+^^^\vl) , (33) 

with the detection flavor state | v^) inEq. (21). 

Neglecting mass effects in the production and detection processes, we approximate 
the production and detection flavor states with the standard ones given in Eq. (2). Then, 
we obtain 

AapiT,L) = Y^KkUpke-'^'^'-'P''^. (34) 

k 

Notice that the consideration of the space-time interval between neutrino production and 
detection allows one to take into account both the differences in energy and momentum 
of the massive neutrinos [22, 23, 24, 25, 26]. 

In oscillation experiments in which the neutrino propagation time T is not measured, 
it is possible to adopt the light-ray T = L approximation (assumption (A3)), since 
neutrinos are ultrarelativistic (the effects of possible deviations from T — L are shown 
to be negligible in Refs. [34, 18]). In this case, the phase in Eq. (34) becomes 

Ek + Pk Ek + Pk ^E 

which leads to the standard oscillation probability in Eq. (8). 

Equation (35) shows that the phases in the flavor transition amplitude are independent 
from the values of the energies and momenta of different massive neutrinos [22, 23, 
24, 25, 26], because of the relativistic dispersion relation in Eq. (4). In particular, 
Eq. (35) shows that the equal-momentum assumption (A2) in section 2, adopted in 
the standard derivation of the neutrino oscillation probability, is not necessary in an 
improved derivation which takes into account both the evolutions in space and in time 
of the neutrino state. 

We have called this derivation of the flavor transition probability "covariant plane- 
wave theory of oscillations" because it is manifestly Lorentz invariant. This is important 
because flavor, which is the quantum number that distinguishes different types of quarks 
and leptons, is a Lorentz-invariant quantity. For example, an electron is seen as an 
electron by any observer, never as a muon. Therefore, the probability of flavor neutrino 
oscillations must be Lorentz invariant [24, 35]. 



FIGURE 1. Space-time diagram representing schematically the propagation of the wave packet of a 
massless particle (a) and the propagation of the wave packets of a superposition of a massless and a 
massive ultrarelativistic particle (b) from a production process P to a detection process D. 

5. WAVE-PACKET TREATMENT 

So far, we have considered massive neutrinos as particles described by plane waves 
with definite energy and momentum. However, the T = L assumption (A3) requires 
a wave packet description. The reason is simple: since plane waves cover all space- 
time in a periodic way they cannot describe the localized events of neutrino production 
and detection. As discussed in introductory books on optics (see [36, 37]) and quantum 
mechanics (see [38, 39]), real localized particles are described by superpositions of plane 
waves known as wave packets. 

Moreover, different massive neutrinos can be produced and detected coherently only 
if the energies and momenta in the production and detection processes have sufficiently 
large uncertainties [40, 41]. The uncertainty of the production process implies that the 
massive neutrinos propagating between production and detection have a momentum 
distribution [21], i.e. they are described by wave packets. 

The propagation of a massless particle between localized production and detection 
processes separated by T ~ L is illustrated schematically in the space-time diagram in 
Fig. la. The interesting case of propagation of a superposition of two neutrinos with 
definite masses, one massless (Vi) and one massive but ultrarelativistic (V2) is illustrated 
schematically in Fig. lb. One can note that in these diagrams both the production and 
detection processes occupy a finite region in space-time, called the coherence region, 
in which the propagating particles are produced or detected coherently. Indeed, the 
uncertainty principle implies that any interaction process / has a space uncertainty ol 
related to the momentum uncertainty by 

a>J~i. (36) 
A point-like process would have an infinite momentum uncertainty and a process with 



definite momentum would be completely delocalized in space. The momentum uncer- 
tainty can be estimated as the quadratic sum of the uncertainties of the momenta of the 
localized particles taking part in the process: 



i^lr-Uv- (37) 

i 

The sum is over the initial particles and the final particles which are localized through 
interaction with the environment. Their momentum uncertainties Op are related to the 

size of their wave packets by uncertainty relations analogous to Eq. (36), 

(38) 

Therefore, the space uncertainty of the process is given by 

It is clear that the particle with larger momentum uncertainty and associated smaller 
space uncertainty gives the dominant contribution. 

The coherence time c/ of an interaction process I is the time over which the wave 
packets of the interacting particles overlap. If the process is the decay of a particle in 
vacuum, the localization of such particle and its decay products is very poor and the 
coherence time oj is of the order of the particle lifetime. On the other hand, if the decay 
occurs in a medium where the decaying particle and its products are well localized or if 
the production process is a scattering process, the coherence time can be estimated by 

i^i)-'-L(-y\ (40) 



where v/ is the velocity of the particle i, because a} must be dominated by the particle 
with smaller ratio o^Vi, which is the first to leave the interaction region. Therefore, in 
general oj > aj, in agreement with the physical expectation that the coherence region 
of a process must be causally connected. 

As illustrated in Fig. 1, one can estimate the size of the wave packet of a massive 
neutrino created in a production process P as the coherence time of of the production 
process, 

(7;~(Tf. (41) 

Let us emphasize that there is a profound difference between the behavior of final 
neutrinos and other particles in the production process. The initial particles have wave 
packets which are determined by their creation process or by previous interactions. The 
initial particles and the final particles which interact with the environment contribute 
to the coherence time of through their contribution to the momentum uncertainty in 
Eq. (37). An initial decaying particle contributes directly to the coherence time of 
with its lifetime. On the other hand, neutrinos are stable and leave the production 
process without interacting with the environment. Therefore, they do not contribute 



to the determination of the coherence time of and the size of their wave packets is 
determined by of. 

Considering now the detection process D, if there is only one particle propagating 
between the production and detection processes, as shown in Fig. la, the coherence size 
of the detection process is determined by Eq. (37), with the sum over all the participating 
particles which interact with the environment and the propagating particle, which is 
described by a wave packet. In the case of neutrino mixing, the neutrino propagating 
between the production and detection processes is in general a superposition of massive 
neutrino wave packets which propagate with different phase velocity, as illustrated in 
Fig. lb. In this case, in the detection process, the wave packets of different massive 
neutrinos are separated by a distance Ax = AvT, where Av is the velocity difference. If 
the source-detector distance is very large, the separation of the massive neutrino wave 
packets at detection may be larger than their size, leading to the lack of overlap [42]. 
In this case, the effective coherence size of the neutrino wave function at the detection 
process is 

<eff ~ ^Ji<yn' + i^r - 7(c7P)2 + (Avr)2 . (42) 

However, Eq. (39) shows that the particle with smaller space uncertainty gives the 
dominant contribution to the coherence size of the detection process. Therefore, if the 
effective coherence size in Eq. (42) of the neutrino wave function is dominated by the 
separation of the wave packets (AvT ^ of) and there is another particle participating in 
the detection process which has much smaller space uncertainty, the different massive 
neutrinos cannot be detected coherently. In this case, there cannot be any interference 
between the different massive neutrino contributions to the detection process and the 
probability of transitions between different flavors reduces to the incoherent transition 
probability 

pmcoherent^£|^^^|2|^^^|2^ (43) 
k 

which does not oscillate as a function of the source-detector distance. On the other hand, 
if all the other particles participating in the detection process have space uncertainties 
which are larger than effective coherence size in Eq. (42) of the neutrino wave function, 
the different massive neutrinos are detected coherently [41], leading to the interference 
of their contributions to the detection process which manifests itself as oscillations of 
the probability of flavor transitions, according to Eq. (8). 

These considerations show that a wave-packet treatment of massive neutrinos is 
important in order to understand the coherence properties of neutrino oscillations. 

6. QUANTUM-MECHANICAL WAVE-PACKET MODEL 

In this section we present a simple one-dimensional quantum-mechanical wave-packet 
model [23, 43] in which the momentum uncertainties of the states which describe the 
produced and detected massive neutrinos are approximated by Gaussian distributions. 
More complete three-dimensional models in which the neutrino momentum uncertain- 
ties are obtained from a quantum field theoretical calculation of the production and de- 
tection processes are discussed in Refs. [44, 45, 46, 47, 48, 21]. 



Neglecting mass effects in the production and detection processes, we describe the 
produced and detected neutrinos in a — > experiment with the wave-packet flavor 
states 

|v^) -L^ak I dPVkiP) Mp)) , \vf) = If/^, / dpwFip) \Vkip)) , (44) 
% J 1, J 



with the Gaussian momentum distributions 



wlip) = 



1 



-exp 



4(aP)2 



¥np)= 



-exp 



(P-Pk) 
4(aD)2 



(45) 

The average momenta pk of the massive neutrinos are determined by the kinematics of 
the production process. They are the same in the detection process because of causality. 
On the other hand, the energy-momentum uncertainties in the production and detection 

processes. Op and o^, may be quite different. 
The flavor transition amplitude is given by 



Aap{T,L) = {vf\e 



-/Er+/PL| P\ 
I ^a/ 



°^ T^Kk^lik / dp exp 

k 



-iEk{p)T + ipL- 



{p-PkY 



(46) 



with the massive neutrino energies 



Ek{p) = \lp^ + ml. 



and the global energy-momentum uncertainty 

1 1 



+ 



(47) 



(48) 



This expression has a correct behavior from the physical point of view, because the 
smaller energy-momentum uncertainty must dominate in the determination of the total 
uncertainty. On the other hand, the global space-time uncertainty cTc = 1 /2<7p is domi- 
nated by the largest of the space-time uncertainties = 1/2(7^ and = 1/2(7^ of the 
production and detection processes: 



D\2 



(49) 



Since in practice the massive neutrino wave packets are always sharply peaked at the 
average momentum (<7p <S E^{pk) /f^kX we can approximate 

Ek{p) ^Ek + Vk{p- Pk) , (50) 
where Ej^ and are, respectively, the average energy and the group velocity given by 

dEkip) 



■ Ek{pk) 



pl + m^ 



1' 



Vk 



dp 



P=Pk 



Pk 
Ek 



(51) 



With this approximation, the integration over dp in Eq. (46) is Gaussian, leading to 



k 



-iEkT + ipkL- 



4(72 



(52) 



Comparing with Eq. (34), one can notice the additional suppression factor for \L — 
VkT\ > Ox due to the wave packets. 

Finally, integrating the space-time dependent oscillation probability P^piT^L) = 
|A(^l3(r,L)p over the unobserved propagation time T, we obtain, for ultrarelativistic 
neutrinos. 



-2ni 



2e2 



rose 
^kj 



with the oscillation and coherence lengths 



rCOh 

^kj - 



rcoh 
^kj 



|Am2.| 




(53) 



(54) 



The coefficient i^, which is the only quantity in Eq. (53) depending on the production 
process, comes from the general ultrarelativistic approximation [23, 49, 50, 51, 34, 18] 



ml 



mt 



" ^ IE 



(55) 



In the limit of negligible wave packet effects, i.e. for L <S L^^*^ and Ox <S L^-^, the 
oscillation probability in the wave packet approach reduces to the standard one in Eq. (8), 
obtained in the plane wave approximation. The additional localization and coherence 
terms 



2e2 



rose 
^kj 



Pkf = exp 



reoh 
^kj 



(56) 



have the following physical meaning [23, 43, 52, 21, 26, 44, 45, 47, 48, 33, 34]. 

The localization term Pj^y suppresses the oscillations due to Am|y if Ox > L^J^. This 
means that in order to measure the interference of the massive neutrino components 
and Vj the production and detection processes must be localized in space-time regions 
much smaller than the oscillation length L'^y. In practice this requirement is satisfied in 
all neutrino oscillation experiments. 

The localization term allows one to distinguish neutrino oscillation experiments from 
experiments on the measurement of neutrino masses. As first shown in Ref. [40], neu- 
trino oscillations are suppressed in experiments which are able to measure, through 
energy-momentum conservation, the mass of the neutrino. Indeed, from the energy- 
momentum dispersion relation in Eq. (4) the uncertainty of the mass determination is 



Smk" = \J{2Ek8Ekf + {2pkdpky 



1^/lEo, 



(57) 




FIGURE 2. Schematic illustration of the separation of two wave packets with different group velocities, 
produced coherently at L = with widths crj determined by the coherence size of the production process. 
The coherence size of the detection process is assumed to be negligible. 



|Am|^ |, the mass 



where the approximation holds for ultrarelativistic neutrinos. If 5mk^ < 
of Vk is measured with an accuracy better than the difference Amlj. In this case the 
neutrino Vj is not produced or detected and the interference of Vk and Vy which would 
generate oscillations does not occur. The localization term automatically suppresses the 
interference of Vk and Vy, because 



-27t' 



I osc 



(58) 



If the condition 



(59) 



which is necessary for unsuppressed interference of and Vj, is satisfied, as usual in 
neutrino oscillation experiments, the localization term can be neglected, leading to the 
flavor transition probability 



-2ni 



rose 
^kj 



kj 



(60) 



which is a function of the distance L, depending on the oscillation and coherence lengths 
inEq. (54). 

In Eq. (60), each k, j term contains, in addition to the standard oscillation phase, the 
coherence term Ffj^, which suppresses the interference of the massive neutrinos and 

Vj for distances larger than the corresponding coherence length, i.e. for L » 1,^°^. This 
suppression is due to the separation of the different massive neutrino wave packets, 
which propagate with different velocities, as illustrated in Figs, lb and 2. When the 
wave packets of V^ and Vj are so much separated that they cannot both overlap with 
the detection process, the massive neutrinos Vk and Vj cannot be absorbed coherently 
[42, 41]. In this case, only one of the two massive neutrinos contributes to the detection 
process and the interference effect which produces the oscillations is absent. However, in 
general, the flavor transition probability does not vanish. For example, if L 3> for all 
k and j, the flavor transition probability reduces to the incoherent transition probability 
in Eq. (43). 



7. CONCLUSIONS 



We have reviewed the standard theory of neutrino oscillations, highlighting the three 
main assumptions: (Al) the definition of the flavor states, (A2) the equal-momentum 
assumption and (A3) the time = distance assumption. 

We have shown that the flavor neutrino state that describes a neutrino produced or 
detected in a charged-current weak interaction process depends on the process under 
consideration. The standard flavor states are correct approximations of these states in 
oscillation experiments, which are not sensitive to the dependence of neutrino interac- 
tions on the different neutrino masses. 

We have presented a covariant plane-wave theory of neutrino oscillations in which 
both the evolutions in space and in time of the neutrino state are taken into account, 
leading to the standard probability of flavor transitions. In this model, no assumption on 
the energies and momenta of the propagating massive neutrinos is needed. Moreover, 
the derivation of the Lorentz-invariant flavor transition probability is manifestly Lorentz 
invariant. 

We have argued that the time = distance assumption derives from the wave-packet 
character of the propagating neutrinos. We have discussed the necessity of a wave- 
packet treatment of neutrino oscillations for the description of the localization of the 
production and detection processes and the coherence of the oscillations. We have also 
presented a simple quantum-mechanical wave-packet model which leads to the standard 
probability of flavor transitions with additional localization and coherence terms which 
have important physical meaning. 

In conclusion, we would like to emphasize that the insight of the founders of the 
theory of neutrino oscillations led them to the correct standard expression for the flavor 
transition probability. Our more modest task has been to clarify the assumptions and to 
try to improve the derivation hoping to elucidate the deep physical nature of neutrino 
oscillations. 
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